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Prefatory. — It is somewhat surprising that so little mathematical work 
should have been done on the subject of epidemics, and, indeed, on the 
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distribution of diseases in general. Not only is the theme of immediate 
importance to humanity, but it is one which is fundamentally connected 
with numbers, while vast masses of statistics have long been awaiting proper 
examination. But, more than this, many and indeed the principal problems 
of epidemiology on which preventive measures largely depend, such as the 
rate of infection, the frequency of outbreaks, and the loss of immunity, can 
scarcely ever be resolved by any other methods than those of analysis. For 
example, infectious diseases may perhaps be classified in three groups : 
(1) diseases such as leprosy, tuberculosis, and (?) cancer, which fluctuate 
comparatively little from month to month, though they may slowly increase 
or decrease in the course of years ; (2) diseases such as measles, scarlatina, 
malaria, and dysentery, which, though constantly present in many countries, 
flare up in epidemics at frequent intervals; and (3) diseases such as plague 
or cholera, which disappear entirely after periods of acute epidemicity. 

To what are these differences due ? Why, indeed, should epidemics occur 
at all, and why should not all infectious diseases belong to the first group 
and always remain at an almost flat rate ? Behind these phenomena there 
must be causes which are of profound importance to mankind and which 
probably can be ascertained only by those principles of careful computation 
which have yielded such brilliant results in astronomy, physics, and 
mechanics. Are the epidemics in the second class of diseases due (1) to a 
sudden and simultaneous increase of infectivity in the causative agents living 
in affected persons ; or (2) to changes of environment which favour their 
dissemination from person to person ; or (3) merely to the increase of suscep- 
tible material in a locality due to the gradual loss of acquired immunity in 
the population there ; or to similar or other causes ? And why should 
diseases of the third class disappear, as they undoubtedly do, and diseases 
of the first class remain so persistently ? — all questions which immediately 
and obviously present themselves for examination. 

The whole subject is capable of study by two distinct methods which are 
used in other branches of science, which are complementary of each other, 
and which should converge towards the same results — the a posteriori and the 
a priori methods. In the former we commence with observed statistics, 
endeavour to fit analytical laws , to them, and so work backwards to the 
underlying cause (as done in much statistical work of the day) ; and in the 
latter we assume a knowledge of the causes, construct our differential 
equations on that supposition, follow up the logical consequences, and finally 
test the calculated results by comparing them with the observed statistics. 

Apparently the first a posteriori work of importance on epidemics was a 
.communication by Dr. William Farr in 1866 (1), in which he maintained, in 
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connection with an epidemic of cattle plague, that the course of the epidemic 
would follow a curve of which the third difference of the logarithm was a 
negative constant. In 1873-74, Dr. G. H. Evans (2) endeavoured to extend 
this law to other epidemics, but ended by expressing disappointment with his 
results. Quite recently, however, Dr. J. Brownlee has continued the work 
in a series of excellent papers (3-8). In the first of these he said that 
Dr. Farr's curve is the normal curve of probability, but, after studying a 
number of epidemics, found that Pearson's Type IV fitted better, and he 
concluded (p. 567) that " an epidemic is an organic phenomenon, the course 
of which seems to depend on the acquisition by an organism of a high grade 
of infectivity at the point where the epidemic starts, this infectivity being 
lost from that period till the end of the epidemic at a rate approaching to the 
terms of a geometrical progression." He admitted, however (p. 517), that 
" other factors, which are not clear, seem to come into play " ; but added, 
" that the epidemic ends because of the lack of susceptible persons has no- 
evidence in its favour, either from the form of the curve or from the facts "' 
(see also pp. 500 and 501). It is obvious from the mere examination of 
many curves of epidemics that they are often remarkably symmetrical bell- 
shaped curves, which, however, frequently tend to fall somewhat more slowly 
than they rose ; and in a later paper (4) Brownlee emphasises this feature of 
symmetry, and says (p. 2) that " the deduction from this phenomenon is direct 
and complete, namely, that the want of persons liable to infection is not the 
cause of the decay of the epidemic. On no law of infection which I have 
been able to devise would such a cause permit of epidemic symmetry. The 
fall must in all such cases be much more rapid than the rise, though, on the 
contrary, when asymmetry is markedly present the opposite holds." In still 
later papers (6, 7) he gives much evidence to show that measles and smallpox 
have distinct periods of recrudescence. 

So far as I can ascertain, a priori researches on epidemiology were first 
commenced by myself since 1899 in connection with malaria. In 1904 I read 
a paper (9) concerned with the random migration of mosquitoes, a subject of 
vital importance in the theory of malaria, which was subsequently dealt with, 
at my suggestion, by Prof. K. Pearson, whose researches were then employed 
by Brownlee (5). Subsequently I endeavoured to find the malaria equations 
bj a priori reasoning (10), and in the second edition of my book (11) extended 
this method to a preliminary general "Theory of Happenings," employing 
chiefly the Finite Calculus (which is useful for malaria). Quite recently 
I published a very brief note (13) recording further advances; and Brownlee 
added another (7), showing how he arrived at some of his results. 

The present paper is a much more advanced and general development of 
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the work commenced in Section 66 (14) of the second edition of my book 
(11, p. 676). It deals merely with the curves which are theoretically 
obtained when we suppose that the infectivity-ratio remains constant or 
proportional to the number of individuals already affected, while, simultane- 
ously, some of these are constantly "reverting," or losing immunity, and 
while both the affected and the unaffected groups are subject to special rates 
of birth, death, immigration, and emigration — as defined in the following 
section. This thesis does not consider progressive decline of the infectivity ; 
yet we shall see in Section VII, that, contrary to what Dr. Brownlee 
supposed, it does often yield curves which are bell-shaped, nearly symme- 
trical, or, when otherwise, decline more slowly than they rose — that is, yields 
curves of the type frequently found in epidemics and does so without 
demanding any other cause for the decline than the exhaustion of susceptible 
material. 

The thesis is stated in terms which are certainly academical, but into 
which it has been resolved only after very careful thought. It has not only 
been dealt with in considerable detail, but some of the more immediate 
consequences have been followed out— because I think that it must first be 
examined and used as a standard before we can proceed to discuss variations 
due to changes in the infectivity. As will be seen, the mathematics presents 
few difficulties ; but, owing to the nature of the subject, I have thought it 
wise to labour the proofs somewhat more than some readers will require. 

Beyond this point I have not endeavoured to go, because the application 
of the equations to numbers of known epidemics (which will ultimately be 
required) can be made only at great length. The paper is therefore purely 
theoretical; but it is one which I think will be needed for future studies. 
In stating the results I do not at all wdsh to contest the conclusions of 
Drs. Farr and Brownlee, which may quite possibly be ultimately found 
capable of being superimposed upon mine. 

Dr. Farr's communication (1) is not republished in his collected works 
(' Vital Statistics/ edited by N". A. Humphreys, London, 1885), but has recently 
been recovered with some difficulty by Dr. Brownlee. It is interesting to 
note that he seems to have attributed the decline of epidemics partly to the 
attenuation of infectivity, but also partly to " the fact that the individuals 
left are less susceptible of attack, either by the constitution or hygienic 
conditions, than those destroyed." The descriptions of his methods furnished 
by himself and by Evans are, however, almost unintelligible ; but it is evident 
from " the calculated series by law " which he gives, that the third difference 
(not the second, as Brownlee said) of the logarithms of the series is a negative 
constant. 
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Some of Brownlee's equations are also difficult to interpret in terms of the 
numbers of individuals affected. I have therefore thought it best to proceed 
in this paper entirely on the basis of my own previous work. 

Part II (which is not yet finished) will contain Sections VIII (Hypo- 
thetical Epidemics), IX (Hypometric Happening), X (Parameter Analysis), 
XI (Variable Happening), and XII (Discussion). 

II. 

The problem before us is as follows. Suppose that we have a population 
of living things numbering P individuals, of whom a number Z are affected 
by something (such as a disease), and the remainder A are not so affected ; 
suppose that a proportion h . dt of the non-affected become affected in every 
element of time dt, and that, conversely, a proportion r . dt of the affected 
become unaffected, that is, revert in every element of time to the non -affected 
group; and, lastly, suppose that both the groups, the affected and the non- 
affected, are subject also to possibly different birth-rates, death-rates, and 
immigration and emigration rates in an element of time ; then what will be 
the number of affected individuals, of new cases, and of the total population 
living at any time t ? 

CD v 

For the solution of this and the subsidiary problems I have ventured to 
suggest the name "Theory of Happenings." It covers many cases which 
occur not only in pathometry but in the analysis of questions connected with 
statistics, demography, public health, the theory of evolution, and even 
commerce, politics, and statesmanship. The name pathometry (pathos, a 
happening) was previously suggested by myself in antithesis to nosometry 
(nosos, a disease) for the quantitative study of parasitic invasions in the 
individual. 

III. 

(i) Let ndt, mdt, idt, edt denote respectively the nativity, mortality, immi- 
gration, and emigration rates of the non- affected part of the population in the 
element of time dt ; and 'Ndt, Mdt, Idt, 'Edt denote the similar rates among 
the affected part. Then, as argued in my previous writings and as will be 
easily seen, the problem before us may be put in the form of the following 
system of differential equations : — 

dP = (n-m + i-e)dt . A + (N-M + I-E)^ . Z, (1) 

dA = (n— m + i— e— h)dt . A + (N" + r)cft . Z, (2) 

dZ = hdt . A + (-M + I-E-r)d* . Z. (3) 

Here dP consists only of the variation-elements n, m, i, e, N", M, I, E, with 
their proper signs, while dA and dZ contain also the happening -element h 
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and the reversion-element r. Obviously dT = dA + dZ, as it should do, the 
elements h and r disappearing in the summation because they denote only 
change of condition and not of being. The number of the non-affected who 
become affected in the element of time is hdt . A, and, of course, this number 
pass over from group A to group Z — as shown in equations 2 and 3. 
Similarly rdt . Z is the number of the affected who revert in the element of 
time and therefore pass over from group Z to group A. The equations are 
not quite symmetrical, since N, which should appear in equation 3, appears in 
equation 2 — because the progeny of the affected group, namely Ndt . Z, will 
generally be born not affected — an important fact which modifies almost all 
the results. 

(ii) The variation-elements n, m, i, e, and N, M, I, E, may sometimes be 
functions of time, especially if the considered events extend over long 
periods ; but it will quite suffice at present to take them as being constants. 
If we have sufficient data regarding them, we can generally calculate them by 
the methods of Section VIII for whatever small unit of time we adopt: 
otherwise they must be conjectured or assumed to begin with, and be then 
ascertained by a comparison of the integrated equations with known facts. 
It is convenient to write 

v EE n—m + i—e, 

V = N-M + I-K (4) 

In some kinds of happening which have no marked effect on the birth-rates, 
death-rates, and immigrations — such for instance as mild maladies, enlistments 
in civil professions, conversions to religious or political parties, entries into 
unlimited societies or trades-unions, etc. — we may have n = E", m = M, i = I, 
# = E, and, consequently, v = Y. In other cases, quite possibly we may still 
have v = Y, though the individual items may be different. Generally in 
marriages T$>n m , in accidents M > m ; in vaccinations M < m ; in total- 
abstainers 1ST > n and M < m ; in military enlistments and in many diseases 
N<n and M>m; while in certain alarming epidemics, especially cholera, 
plague, and malaria, we find in addition that I < i and E > e — so that in these 
cases, which are particularly our present subject of study, we generally have 
v>Y. In fatal accidents M = 1 and N, I, and E all = 0. When considering 
happenings among the same individuals, we also put N", I and E at zero, 
though M may be anything ; and in other special cases we omit various 
elements. If the surrounding population is not affected, 1 = 0; and if the 
affected individuals cannot move, E = : so that it will be seen that the 
equations can be made to cover a wide field. The theory will sometimes also 
apply to inanimate objects — as for instance in commerce, where the variation- 
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elements may be taken to mean manufacture, waste, exportation and 
importation, though the original equations may have to be somewhat 
modified. 

(iii) The reversion-element may also at present be taken as a constant and 
be calculated for the same unit of time as is used for the other elements. In 
the case of independent happenings (to be defined presently) rclt means merely 
the proportion of affected individuals which become, in element of time, 
capable of being re-affected — as by divorce in marriage, and by discharge in 
certain employments. In dependent happenings, however, it implies also the 
loss of capacity for affecting others. Thus in infectious diseases (which are 
dependent happenings) it implies loss both of immunity and of infectivity — 
not recovery from sickness, which is merely an episode of affectedness from 
disease. In some diseases, such probably as leprosy and organic diseases, 
r is zero or nearly zero ; in others with long-continuecl immunity, such as 
many zymotic diseases, it is low ; and in others again, with comparatively 
quick loss of immunity and infectivity, such possibly as nasal catarrh or 
dengue, it must be larger. In many diseases, however, it is quite unknown 
at present and must in fact be calculated from the integrated equations 
(which, it is hoped, will prove of use for this very purpose). In the case of 
slight accidents it is unity ; in fatal accidents it is zero ; and in snake-bite or 
heat-stroke, where recovery is quick if it occurs, it is high, though at the same 
time the death-rate M is also high. In controversial parties due to rational 
divergence of opinion it should be high ; in ordinary party-politics it is, in 
fact, very low. 

(iv) The most important element is the happening-element, h. We should 
clearly understand that in most cases the happening, whatever it is, does not 
select only the non-affectecl, but tends to fall on both groups alike. If, 
however, it chances to fall upon individuals who are already affected, 
it merely re-affects them and does not cause them to pass from, one group to 
another. Eeally the total number of individuals to whom the happening 
occurs in element of time is licit . P, that is, licit . (A + Z). But the individuals 
numbering lidt.Z do not count because they are already affected, and there- 
fore do not appear in the equations. On the other hand, the actual number 
of new cases (which we may denote by Ydt) is hclt . A ; and this quantity 
does appear in the equations 2 and 3. Or we may obtain F independently 
from the proportion 

mi a 



hclt . P P' 

so that, again, F = hA = h (P — Z). This is an important subf unction in all 

cases. 
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(v) Different kinds of happenings may be separated into two classes, 
namely (a) those in which the frequency of the happening is independent of 
the number of individuals already affected ; and (6) those in which the 
frequency of the happening depends upon this quantity. To class (a) belong 
such happenings as many kinds of accidents and non-infectious diseases due 
to causes which operate, so to speak, from outside ; and to class (6) belong 
infectious diseases, membership of societies and sects with propagandas, 
trade-unions, political parties, etc., due to propagation from within, that is, 
from individual to individual. In the former case, h or F will be constants^, 
in the latter case they will be functions of Z. 

To proceed now to the integration of equations 1, 2, and 3 — we may 
observe that the solutions, though they belong to the same class of functions, 
differ specifically in different cases. The integrals are easy to obtain by 
ordinary methods* ; but the most elegant and useful method is to put all 
the cases into similar forms which can be brought by substitutions into the 
same immediately-integrable differential equation of the simplest type. We 
will first take a rapid survey of each solution as it is obtained. 

IV. 

(i) Independent Happenings. — Here h or F is supposed to be constant. 
If h is constant, the happening falls on the same proportion, hdt,' of the 
population in every element of time. 

Putting P — Z for A, using v and V for the variation-elements as in 4, 
and setting x = ZP, equations 1 and 3 become 

dP/dt = vT-(v-~Y)xP, (6) 

dxT?/dt = kP(l-x) + (Y--N-r)xP, (7) 

dx~P/dt = xdP/dt -f Tdx/dt. 

Eliminating dxP/dt and dP/dt from these, we find that P cancels out also, 
so that we have 

dx/dt = h-(h + v-V + N + r)v + (v-~Y)x 2 . (8) 

This has one integral form if v = V, and another if v •=£ V. 

(ii) The Equivariant Case. — Here the happening is such that it does not 
affect the sum of the variation-elements of the affected group — though, as 
already mentioned, the separate variation-elements may be different ; that 
is, the total population is not altered. An example is conversion to some 
philosophical creed without propagandism. Another example is the 

* The equations are familiar in connection with many statistical, chemical, and physical 
problems. 



212 Lieut. -Colonel Sir R. Boss. 

happening of slight accidents — in which, however, r is generally unity ; and 1 
a third example is the happening of a certain standard of wealth, which, we- 
assume, tends to diminish the birth-rate, death-rate, and immigration and 
emigration of the affected, each by an equal decrement, simultaneously. In 
this last case, however, as the children of the affected will tend to be born 
affected, we should remove N" from equation 2 to equation 3, and omit it 
from equations 7 and 8. In all these cases v = V. 

(iii) To integrate equation 8, however, write it in the form 

dx/dt = K(L-x), (9> 

where K = h + N + r and L == h/K. Observe that the tangential, dx/dt y 

vanishes when x = L and becomes negative when #>.L. Put y =E L— x\ so- 

that x = L— y. Then 

dx/dt = —dy/dt = Ky, 

1/y . dy = -Kdt. (10) 

Integrating both sides of this, 

log y = — -K£ + constant. 
If y Q is the value of y at the beginning of the happening, 

log y Q = constant, 
and y = yo€~ Kt . (11) 

Therefore (L-— x) = (L— Xo)e~ Kt , 

and ^ = L~-(L-^o)€~ K *. (12> 

This gives the proportion, %, of the total population who are affected at the* 
time t, this proportion being x when t = 0. 

To find P, put v— V = in equation 6, and we obtain at once a differential 
equation of the same form as equation 9 : so that 

P = P e rt , (13> 

where P is the total population at the beginning of the happening, when t = 0. 
This is an important function which expresses the natural increase of the 
whole population due only to the natural variation-elements, v. Of course v 
is small when the increase of population is not very rapid. When v = 0, the 
population remains constant, the births, deaths, immigrations and emigrations 
annulling each other. 

To find Z we have Z = xP. 

In III (iv), we denned F to be the number of new cases and showed that. 
F = h (P - Z). If / = F/P, we have, when x = 0, 

/= h(l-x) = A(1~L + Le- K *). (14) 
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The functions for Z and x give what may be called respectively the actual 
and proportional curves of affected individuals; and those for F and/ give 
respectively the actual and proportional curves of new cases. 

To find t for any assigned value of x, we have 

Kt = log e (L-^ )~log € (L-^) == 2-3025851 log 10 [(L-^ )/(L-^)]. (15) 
The following equations are required for analysing the curves x and/: — 

cPx/dfi = -K 2 (L-aO> (16) 

d//tf* = -AK(L-a?)(l-a?), (17) 

d»//d^ = AK 2 (L - x)(L + 1 - 2a?), (18) 

(iv) In order to analyse the curve a; (of affected individuals) we first observe 
that as h, N", r are always positive, K and L are also always positive and L is 
always less than unity, while e~ Kt diminishes and finally vanishes as t increases. 
Hence when t is very large, x reaches the limit L (equation 12) and never 
exceeds it. Its tangential (equation 9) begins at the value KL( = A) when 
Xo = 0, always remains positive, and gradually diminishes to zero when x = L 
and t = oo . Its second tangential (equation 16) is always negative. Hence 
as t increases from zero to infinity, x also increases, but with a diminishing 
increment, until it approximates to the asymptote L ; and it is always convex 
to the axis of t. 

On the other hand the curve / (equation 14, of new cases) begins at its 
greatest value h (when x = 0) and constantly diminishes, its tangential 
(equation 17) being always negative. Its second tangential is always 
positive, since J(L+1) is greater than L ; and the curve is, therefore, convex 
to the axis of t, and, when t is large, approaches its final value h(l — L), that 
is (N-hr)L, which is always less than h. Let this constant be denoted by I. 

It may have been thought that, at least in the case of constant happenings, 
the whole population would become affected; but it has been shown that 
this is not the case, because x cannot exceed L, which is a proper fraction. 
What is it then that limits the value of x ? Equation 8 may be written 
(when v = V) 

dx/dt = h(l - x) - (N + r)x = /- (N + r)x, (19) 

and dx/dt vanishes when/= (N + r)x. When it vanishes, x = L and/= I; 
that is, the proportion of affected individuals comes to a limit when the 
proportion of new cases exactly balances the proportion of recoveries and 
births (it being assumed by supposition that the progeny of the affected are 
born not affected). If, however, N" and r are zero (and only then), L = 1, 
and the whole population does become affected in time. It should also be 
observed that I . = A(l— L) =&Al/P, where Al is the final number of 
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non-affected individuals when x = L and £ is very large ; that is to say, even 
after the limit of x has been nearly approached, the happening continues and 
constantly affects the same proportion of the non-affected population as 
before — but now A also comes to a limit, namely P(l— L). 

The easiest way to assess the increase of x as a function of t is to give to 
e Kt the successive values 1, 2, 3, 4 ... 10, 100, 1000 .... Then it is clear from 
equation 12 that, if x = 0, x will have the corresponding values 0, |L, fL, 
|L, ... 0'9L, 0*99L, 0-999L, and t will have the values 

i*!±^ log 10 (l, 2, 3, ... 10, 100, 1000 ...). 

An important point is reached when e Kt = 10, for then 

£ = ~ 2-302... = T (say). (20) 

At this point x has increased to 9/10ths of its ultimate value L — that is, 
most of the change of w 7 hich it is capable has been effected, and it increases 
only slightly afterwards, namely only by 0*09 and 0*009 when t = 2t, and 
?= 3t, and so on. On the other hand, x reaches half its ultimate value very 
quickly when t = ^r nearly. 

However large t may be, e~ K * never quite vanishes, and therefore x never 
quite reaches L. It is convenient, then, to find the value of t when the 
number of affected individuals reaches the limit less one individual — that is, 
when Z = LP — 1. This is obtained at oncer from equation 15, since 
'L/(L-—#) ='LP/(LP — Z); and the required value is 

t = t logioLP = t (say). (21) 

For example, if L = J and P = 2,000,000 and remains nearly constant, the 
number of affected individuals is 900,000 when t = r, and is 999,999 when 
t = 6t. 

The ratio xji is the tangent of the radius-vector from the origin to the 
point (x, t). If xq = the value of this tangent is h, that is, dx^jdi, when 
t = ; and is 9L/10t, that is, 0*39086/^, when t = r ; so that, as the curve of 
x increases from zero to nine-tenths of its ultimate value, it lies wholly 
within the angle formed by these vectors. And the magnitude of this 
angle depends only upon k, and not upon the other elements, 1ST and r. 

The constants L, I, and t may be written in detail, 

t 1 j J- Z'ovZ. . . /99\ 

"l + (N + r)A- lf u ""7r 1 + (N + ? i )- 1, T "A + N + r" K " } 

From these equations it will be easily seen that if h increases while 1ST and r 
remain constant, L and I will increase and t will diminish ; but if h remains 
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constant while KT or r or both increase, L and r will diminish and I will 
increase— just as may be expected. If r = 1 and N^O, I > L, since 
I = (N-fr)L; that is, as all the cases recover as soon as they occur, the total 
number of cases will not exceed the new cases during unit of time — as 
happens with slight accidents such as mosquito bites, where there is nothing 
in one happening to prevent a second or third occurrence of it. If 
N-fr = 0, L = 1, I = 0, and r = h, and if h = IsT + r, L = J, I = \h, and 
r = 1*151.... 

(v) It remains to examine some questions which require for their solution 
the integrals of P, Z, A, and F — all easily obtained. 

From equation 13 we have at once 

\Vclt = lT 6 vt dt = P/-U + constant, 

Ppd* = (P-P )/t>. (23} 

To find the integrals of Z, A, and F, we integrate both sides of equation 7 

as there written, or as written in the form dz/dt = &P — (K — V)Z, and 

obtain 

Z = \Edt + ( V - N" - *■) Jz<ft, (24) 

or = AjP<«-(K- V)JZd* ; ( 25 ) 

whence Zd* = J =— ^ - J ^-~, (26) 

Jo v h 



1 — io j Li — Zjq 

V k 



Mt = (1 - J) i — - + J £i - r £ ^ (27) 



where J = h/(K~Y) and 1 — J = — (V— N"-fr). And, of course, 



! 



F^ = h Adt 
o Jo 



The same equations can be readily obtained by direct integration from 
equations 12, 13, etc. 

The concrete interpretation of these definite integrals of P, Z, and A, is 
that they express the total number of time-units lived during the period t by 
the whole population, by the affected population, and by the non-affected 
population respectively. 

For example, suppose that a population P has doubled itself in 1000 
time-units, then, since 1000^ = log e 2, v = 0*00069315, and the total number 
of time-units lived by all the individuals together is Po(2 — 1)1 v — that is, 
1442*7 Po time-units. 

It should be observed that if v is very small, equation 23 takes an 
indeterminate form, the value of which is the value of P (€ vt — l)/v when 
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v = — namely, £P . That is, when the population remains nearly constant, 
the total number of time-units lived during t is £P , as we should expect. 

The same integrals divided by t will give the mean numbers of the total, 
the affected, and the non-affected populations respectively during the period t. 
Thus, if the happening is such that it does not affect the normal rate of 
increase or decrease of population, then the mean population will be 
(P — Po)/log e P/Po — as we should expect. 

When multiplied by the appropriate constants, these integrals will also 
give the total number of the variation-events, births, deaths, immigrations and 

emigrations. Thus the total progeny of the non-affected will number n Adt, 

Jo 

and the total deaths among the affected will number M Zdt. The sum of 

Jo 

the natural variation-events in the whole population is (%— m + i — e)\ T?dt, 

Jo 

that is, P— P . Similarly r Zdt gives the total reversions : and the integrals 

Jo 

multiplied by h give the total number of time-units lived in which the 
happenings have occurred — which is the same as the number of individuals 



to which the happenings have occurred. Thus h 



't 



Adt is the total number of 

o 

happenings (not necessarily for the first time) among the non-affected — that 

rt 



is the total number of new cases, 



Fdt. 

Jo 



It is now easy to interpret equations 24, 25., 26,-27. Since V— N— r 
= 1— E— M— r, equation 24 means that Z, the total number of affected 
individuals remaining alive at the time t, is equal to the sum of all the new 
cases plus the affected immigrations and less the affected emigrations, deaths, 
and recoveries. Also as K— V = A + r + M + E— I, equation 25 means that 
the same quantity Z is equal to the total happenings plus the affected 
immigrations, less the affected emigrations, deaths, and recoveries, and also 
less the happenings among the affected— -which, as stated in III (iv), do not 
count. 

V. 

(i) Independent Happenings : the General Case : v ^ V. — The method of 
working employed for the equivariant case has been given at some length 
because it will serve also for the other cases. We now proceed to integrate 
the general equation 8, namely 

dx/dt = h—(k+v—. Y+N + r)x + (v— Y)x 2 . 
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Write this in the form 

dx/dt = K(L-x)(L' -x), (28) 

where K = v—V and L and 1/ are the two roots of the quadratic in x. Let 
L = *-/3 and 1/ = * + /3, where a = (h + "K + r + K)/2K and /3 2 = ct 2 -h/K. 
These roots are always real because 

(& + N + r + K) 2 -4AK = (h-K) 2 + 2(h + K}(N + r) + ($ + r) 2 , 

of which the left side is positive when K is negative, and the right side is 
positive when K is positive. And the roots are both positive when v > V. 
Observe that dx/dt vanishes when x = L or I/, if this ever occurs, 
(ii) Put 

y = (I/-x)J(L-x), a? == L + (I/-L)/(y-l). (29) 

dy/dt = K(L'-L)y, y = yotpw-v*. (30) 

as in equation 11. Therefore, as I/— L = 2/3, 



or 03 = L — 



L— a; L— a3 
2/8(L-ar ) 



(31) 



(L'-aJbJe^^-CL-aro) 

= L { 1 -. LJ^-L } Wheaa;o = () - (32) 

The first tangential of this is given in equation 28 ; the second tangential 

is 

d 2 x/dt 2 = -2K 2 (x-x)(L-x)(L'-x). (33) 

(iii) The analysis of x, the proportional curve of affected individuals, depends 
upon whether K, that is v— V, is positive or negative. If it is positive, we 
are generally concerned with cases of Injurious Happenings; because, as 
explained in Section III (ii), the nativity, mortality, immigration and 
emigration are likely to be respectively less, greater, less, and greater among 
the affected in such cases- than among the non-affected. Non-infective 
diseases should be examples in point. On the other hand, if v—Y is 
negative, we shall be generally concerned with cases of Beneficial Happen- 
ings, which improve the natural variation elements of the affected — such for 
example as conversions to total abstinence are said to be (if they do not 
diminish i and increase e in places where they prevail, such as the Prohibition 
States in America.) 

If K is positive, e 2K ^ in equations 31 and 32 increases indefinitely with t 
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(since yS is the positive value of the radicle in the roots L and I/) ; so that 
x approaches the limit L, that is, a — /3 ; and never exceeds it, and therefore 
never reaches the greater root I/, that is, a + /3. Therefore dx/dt is always 
positive for values of x between zero and L ; and d 2 x/dt 2 is always negative 
for the same values, since a > a — j3. Hence in this case x has a form similar 
to what it has in the equi variant case. That is, it begins (when x = 0) at 
zero at an angle of which the tangent is h; constantly increases with a 
decreasing increment, so that it always remains concave to the axis of t; 
and finally approaches the limit L when t is very large. 

When K is negative, the quadratic expression in equations 8 and 28 has 
only one change of sign, and therefore only one positive root, I/, or <* + /3. 
The other root, L, does not concern us because we consider only positive 
values of x. In this case, e 2K ^ diminishes indefinitely as t increases, and 
x finally reaches the limit I/. Here, too, dx/dt always remains positive for 
possible values of x, that is, while t varies from zero to positive infinity ; 
so that x again always increases from zero to I/. If h + v— V-f N" + r, that 
is, 2K&, is positive, d 2 x/dt 2 is always negative in this case also, and x is 
therefore concave to the axis of x as before. But if this expression is 
negative, d 2 x/dt 2 will be positive at first, but will change sign when x = <*, 
and finally vanish when #=1/ ; so that, in this case, x will be first convex 
and then concave to the axis of t. This is the only difference in general 
form. It is seldom that v— V can be greater than fo + N+r, which is 
always positive ; hence x has the same form in most cases where h is a 
constant. 

(iv) For the proportional curve of new cases,/, we have 

/=A(l-a); (34) 

df/dt = —h dx/dt d 2 f/dt 2 = -hd 2 x/dt 2 (35) 

Hence, when x = 0, / begins at the value h and, as t increases, gradually 
falls to the limit A(l— L) when K is positive, and to the limit A (1—1/) 
when K is negative. The curvature is the opposite of that of x- — as in the 
equivariant case. 

To find the value of P, the total population, write equation 6 in the form, 

1/P • 02/ dt = v-Kx, 
Therefore log P = vt— \Kxdt 

J m— e~ 2K ^ 
= (v — KL)£ + log(y Q — e~ 2K ^) -f constant. 
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On evaluating the constant we find that we can arrange P in either of two 
ways, namely, 

P = p o6 ^ gg 6 6 -2KL^ Qr p == p o6 rt ,ygg P . -V ^ (36) 

yo— 1 3(0—1 

If xo = and £ is large, these become, when K > and K < 0, respectively, 

P = p oe ^ *±_ e -KLe or p = p oe * „z^ e ~~km (37) 

The former expressions are suitable When K is positive, and the latter 
when it is negative. Thus in the case of injurious happenings when K 
is positive, the population will diminish indefinitely unless the natural 
increase denoted by v is large enough to compensate for the decrease, KL, 
due to the happening (L being positive in this case). In the case of 
beneficial happenings when K is negative, L is also negative while L' 
remains positive, so that the happening enhances the natural increase of the 
population due to the happening. 

, As with equivariant happenings, x never quite reaches either L or 1/ ; 
and it will therefore be useful to find the value of i when the effected 
population equals the limit less one individual, that is, when Z = LP— 1 
or I/P — 1, or when x = L— 1/P, or L'-~l/P. This can easily be ascertained 
from equations 29 and 30, and we have to find t from the equation 
1/P = (L' — L) /(?/—■ 1) when K is positive, and from -> the ! equation 
1/P =r y(I/ — L)/(y — 1) when K is negative. From these we have when 
Xo = 

(KL'-v)t = log € LP (K>0) 

(KL-v)t = log.LT . (K<0) (38) 

To find the integral of P, we write P in the form (when x -= 0) 

P = _ * ° P (»-KL)I l>Po Jv-Kb^L 

Ju — i_j J_j — Xj 



Therefore 

•* 1/ p -KLt__ p -Kh T P -KI^__,,-KL' 

X.LW — • lot YT t tFT — ' - 1 € 7";" ■ v""~" Ty Y , — '. {ovl 



i: 



The integrals of Z and F are found from this by the same methods as were 
used in equations 24 and 25, or by integrating both sides of equation 6. 
From the latter we have 

(v-Y) [z* = v [Pd»-(P-Pe). (40) 

The remarks made in Section IV (v) apply here also. 

VOL. xoil— A, s 
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VI. 

Independent Happenings : Constant New Cases.— A. third example of inde- 
pendent happenings remains to be considered, namely, the one referred to at 
the beginning of Section IV, in which F, the number of new cases, is always 
a constant. The only instances of this which I can recall are the cases of 
certain societies (such as the Eoyal Society) which elect a fixed number 
of new members every year. Indeed, this process of increment seems to be 
too artificial to be seen in nature, unless, perhaps, in certain cases of seeding 
or spore formation, and we need only note the form of the principal function 
for future reference if required. 

Here, since F is a constant, let F = c. Then equations 6 and 7 

become 

dP/dt = *>P-(t;-V)Z, (41) 

dZ/dt = c + (V-N-r)Z, (42) 

= K(L-Z). 

where K =E M — I + E + r and L = c/K. The solution of this is given in 
equation 12, except that Z is here substituted for x. If Z = 0, we have 



Z = c 



M — I + E-f r ' 



It- is to be hoped that the case of the Eoyal Society is at least an equi- 
variant one. 

VII. 

(i) Dependent Happenings: Proportional Happening. — Eef erring to the 
definition of dependent happenings in Section III (v), we see that in such 
h must be a function of Z, and consequently of t. First consider the case in 
which each affected individual affects or infects c other individuals in unit of 
time, c being a constant. This may be taken as being a first approximation 
to the study of those very important happenings, the infectious diseases. 

The total number of the happenings which occur in element of time will 
be eZdt\ but, as in the preceding cases, some of these may chance upon 
individuals who are already affected, and the number of new cases in element 
of time, namely, F dt, will therefore most probably be given by the proportion 

Ydt/eZdt = A/P, (43) 

that is F = cZ(l— x), and h = ex. 
(ii) Equations 1 and 3 now become 

dF/dt = v P - (i? - V) a?P, (44) 

(&P/^ = cajP(l-ic) + (V-N"-r)icP. (45) 
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Treating these as in Section IV (i), we find that, as before, on eliminating 
dxP/dt and dVjdt, P also cancels out, and we have 

dx/dt = (c— v + V — N" — r)x— (c— v + Y)x 2 > (46) 

= Kx(L-x), (47) 

where K = c—v + Y and L = 1 — (lST + r)/K. Here again d#/dft vanishes 
when # = L. Integrating by substitution as before, put 

y = (L—x)/x f # = L/(l+y). (48) 
Hence - ^ -^ ^ = K ,— ( L - -^ = KL y 



(1+y) 3 * l + y\ 1 + 2// (1 + 2/) 2 

rfy/y = -KL* 2/ = ijoe- KLt (49) 



L 

re = 



1-f (L/x ~l)e~~ KLt 

JjXq€ ku 



(50) 



X € KU + (L — Xq) 

Before analysing this it is advisable to obtain the integrated expressions 

for P and Z. Dividing equation (44) by P and integrating both sides, 

we have 

logP = vt—(v — Y)\xdt 

= vt—(v—Y)/K . log (# e KL * + L — xo) + constant, 

L \ (»-V)/K 



whence P = P e^ ( ,^ — - ) (51) 



and Z - Z ^ KLt {^4l^r ' (52) 

where the fraction within the brackets obviously = (x/x ) e~~ KLt . 
The explicit value of t is given by 

KL* = log e ^-log e ^. (53) 

(iii) For analysing the proportional curve of affected individuals, x, the first 
tangential is given in equation 47. The second and third tangentials are 

cPxfdt 2 = K 2 x(L-x)(L-2x), (54) 

d 3 x/dt s = K%(L—^)(L 2 — 6L^ + 6# 2 ). (55) 

From these and from equations 50, it will be seen that x and all its 
tangentials are zero when x = and = L. If x Q is very small, the curve 
begins at nearly zero when t = ; and then increases very slowly at first 
and more rapidly afterwards, remaining convex to the axis of t to begin with. 

s 2 
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When, however, x = JL, there is a change of curvature and x becomes 
concave to the axis of t, and then finally approaches the ultimate limit L 
when t is very large, and never exceeds that limit. Its maximum rate of 
increase, namely |KL 2 , is reached when x = |L ; and then 

KU = log[(L—^o)/^o]. 

The curve is now seen to be a symmetrical one with this point as the centre 
of symmetry. Moving the origin to the point by putting x = ^L-\-x f and 
t = 1/KL. log[(L — ^ )/^o] + ^ 5 we obtain 

x' = JL(l~e- KL 0/(l + 6~ KLt ). (56) 

Substituting — t' for £' in this, we find that x' also merely changes its sign 
without changing its numerical value. In fact, x has the general shape of a 
long-drawn-out letter S. 

(iv) The curve dxjdt is also a symmetrical one with its centre of symmetry 
at the same distance from the primary origin, namely log[(L—^ )/^o]- When 
Xo is small, the curve increases slowly at first, then rapidly ; reaches a 
maximum value, |-KL 2 , at the centre of symmetry, and then falls just as it 
rose. It has two changes of curvature which are the two roots of the 
quadratic expression in d 3 x/dt s , namely when x = (i + v^T^^J ^hat i s wnen 



x 1 = + 0-28867L, 



and KW = + log [(^3 - l)/( v / 3 + 1)] = + 1-31697. (57) 

At these points the ordinate of dxjdt, that is, of K.x(L~-x), is 

KL (| + v tV)(^~~2"~*v tV)> 

that is f -|-KL 2 , or two-thirds of the maximum ordinate at the centre of 
symmetry ; and the ordinate is, of course, the same at both points of change 
of its curvature, as also at other points equidistant from the centre of 
symmetry. 

The curve of dxjdt is therefore a regular bell-shaped curve, very similar to 
those often, found in epidemics. 

By the conditions of this kind of happening, x and Z can never be zero, 
for if they were there would be no new cases. It is convenient therefore to 
take Z as being one individual— that is, Z = 1, x = 1/Po, and dxjdt = 
K(P -l)/Po 2 , or nearly K/P . 

As with the kinds of happening previously considered (see equation 21), 
x and Z never quite reach the limits L and LP respectively ; it is useful then 
to find when the number of affected individuals reaches the limit less one 
individual, that is when Z = LP— 1 and £ = L — 1/P. If Z = 1 and 
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% = 1/P as just suggested, then, owing to the symmetry of the curves, this 
figure will be reached when x = L— Xq ; that is (equation 53) when 

Kli = log[(L— xq)/x ] — log[^ /(L— xq)], 

= 2 log (LP - 1) = 2KLt (say). (58) 

Thus if the total population is large and L is not very small, r, the 
abscissa of the centre of symmetry when x = 1/Po, is very nearly equal to 
1/KL.logLPo; and at twice this period from the beginning the limit of 
affected individuals less one will be reached. We observe also that when 
t = 2r the ordinate of dxjdt will be the same as it was at the beginning, 
namely nearly K/P . 

If xq = 1/Po, the tangent of the radius-vector from the origin to the 
summit of dxjdt is JKL 2 /t, that is 

K*L»/41og(LP -l). (59) 

From equations 48 and 49 we have 



e 



KLt 



xJj—Xq xP LPq — A'oPo Z LP — Zq 



(60) 



a?o L—x . xqPq LP— ^P Z LP — Z 

This provides a useful formula for a cursory estimation of the increase of Z 
with respect to t when the latter is small. For if P and P, Z and Z are 
not very different, and LP is large and Z small — as at the commencement 
of the happening — then this equation is nearly 

Z = Z e KLt . 

Suppose that Z is one individual (t = 0), then Z = 2, 3, 4, . . . individuals 
as € KU increases through the values 2, 3, 4, ... ; that is, as t increases through 
the values 1/KL . logio (2, 3,4, ...). This same thing does not hold true, 
however, when Z is larger. 

Similarly, we may have at first F = cZ e KL *— — c(Zq€ ku ) 2 . 

1 o 

(v) The actual and proportional curves of new cases are respectively: — 

F = cZ(l — x) t 

f — cx(l—x). (61) 

To analyse / we have 

df/dt = cTte(L-x)(l-2x). (62) 

d 2 f/dt 2 = cK a ^(L-a?){L-2(2L + l)x+6x 2 }. (63) 

And, since L = 1— (N+r)/K, we may also write /in the form 

Kf/c = dxjdt + (N + r)a. (64) 

The form of the curve / depends upon the value of L. If L = 1, that is, 
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if (N + r)/K is zero or very small, the curve /equals or approximates to the 
curve dxjclt multiplied by the constant c/K — which is generally greater than 
unity. That is, in this case / is a symmetrical bell-shaped curve with the 
properties just described. In any case, however, / is such a curve plus the 
constantly-increasing function c(N + r)x/K, or c(l — li)x. We may, therefore, 
expect that, at first, when x is small/ will follow the graph of dx/dt, but that 
as x increases, the ordinates of /will become increasingly greater than those 
of dx/dt. At last, when t is very large and x approximates to L, dxjdt 
becomes nearly zero and / approaches the constant value cL(l — L), at which 
it remains indefinitely. Thus / is quite symmetrical only if L = 1 ; but is 
nearly symmetrical if L is only a little less than unity. 

It is evident from equation 62 that clfjdt vanishes when x equals either 
L or J. That is, if L > J, / reaches a maximum when x = |-, the ordinate of 
this maximum being exactly c|(l — J), or |c; and as x increases above this 
value, / then decreases, and falls towards its ultimate constant value 
cL(l-— L) as x approaches L. If L = \ exactly, df/dt vanishes only once 
(and that not quite), namely, when x = L = J, so that / now always increases 
until it reaches its maximum, \c — which is also its ultimate value. If 
L<^, however, x never attains the value J, and consequently df/dt has no 
vanishing point for any values of x between Xq and L, and therefore again 
always increases for all considered values of x and therefore of t. Hence the 
curve of /has two forms: (Type I) if L>|, an irregular bell-shaped form 
(becoming regular if L = 1) ; and (Type II) if L>-|, a drawn-out S-shaped 
form somewhat similar to that of x itself, but not usually symmetrical. 

From equation 63 it will be seen that / may have two changes of 
curvature as x varies from x to L, namely, at the two roots of the quadratic 
expression in the value of d 2 x/dt 2 . The roots are 

6« = (2L+l)±v / [(2L + l) 2 -6L]; (65) 

and are real and positive. The lesser root is always less than L, and is 
therefore always attained by x as it varies from x Q to L — so that / always 
has at least one change of curvature. The greater root is also less than L 
when L>-| ; is equal to L when L = -| ; but is greater than L when L< | — 
so that the second change of curvature occurs only when L> J. (As x never 
quite reaches L, there is no second change of curvature when L = |.) These 
results are therefore as were to be expected from the two forms of the curve 
off already discussed. 

As already seen, when L>|, the summit of the /-curve is \c and is 
reached when x = |. The corresponding value of t is 

t = T-l/KL.log(2L-l), (66) 
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the logarithm being negative since 2L— 1, that is, l-~2(N-f r)/K, is less 
than unity. 

It has also been seen that when L > |-, / falls after reaching its maximum, 
\c 9 to a constant value, cL(l-L). The difference between the maximum 
and the ultimate values of / is therefore 

and their ratio is 1/4L( — L) = X (say). (67) 

As with independent happenings (e.g. equation 22), we shall use the 
symbol I to denote the limit which / approaches when t is very large, and 
this will be examined further in Section VIII (iii) (Part II). 

(vi) The curves of new cases in proportional happenings are especially 
important because, rf the original assumptions upon which they are based are 
sufficient in themselves to explain time-to-time variations in the frequency 
of infectious diseases, these curves should agree, at least for a first approxima- 
tion to the truth, with the curves of new cases actually observed in such 
diseases. Now, according to Brownlee (Section I), the curves of many 
epidemics, especially of the short and sharp zymotic diseases, tend to be 
remarkably symmetrical hell-shaped curves, roughly similar, in fact, to the 
/-curve developed from certain values of the constants. 

We now proceed, therefore, to examine these constants with greater care. 
They are in detail, 

K = c-t> + V, KL = c-h V-N-r, L = l-(N + r)/(c-0 + V) (68) 

The particular elements have been already discussed in Sections III and V. 
The natural-variation element v governs the natural increase of the popula- 
tion apart from exceptional happenings ; but V, N", r are connected with the 
special happening under consideration. If the period covered by such 
happening is short compared with the average life of the individuals con- 
cerned, v may be so small as to be negligible, and N" also may be very small. 
In injurious happenings such as infectious diseases, v— V should be positive 
(Section V (iii)), and V will be negative if M>N" and E>I, so that v—Y 
still remains positive even when v = ; but, as a rule, except in very fatal 
diseases, v— V is very small. The element r may be unity when reversion is 
immediate, as in slight accidents, but may be nearly zero when the immunity 
conferred by one attack of a disease lasts for a long time on the average, as 
in many zymotic diseases. Both N" and r are only positive quantities, so 
that v— V + N + r is usually a small or very small positive fraction in the 
case of infectious diseases (its actual magnitude depending, of course, on the 
magnitude of the time-unit taken as well as on the elements themselves). 
The element of special importance is c, that is the average number of 
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individuals infected or reinfected in unit of time by each already infected 
individual. But while v, V, N", r may be taken as fixed for the kind of 
happening considered,; c is generally unknown, and, indeed, one of the 
ultimate objects of such studies as these will be to ascertain its value. At 
present, however, we take it to be a constant. 

We observe first from equation 53 that if KL is negatives must constantly 
diminish as £ increases, which means that a change must have occurred in at 
least one of the original constants. Hence we assume at present that KL is 
is positive, and reserve the study of the case when it is negative for 
Section IX, (Hypometric Happenings; Part II). 

Secondly, as KL is taken to be positive, c, whatever it is, must be greater 
than v— Y 4- N-f-r— generally a small positive fraction ; and therefore greater 
than each item individually. 

Thirdly, as ov-V + N + r and 1ST and r are positive, c> v — Y ; that is, K is 
always positive. 

Fourthly, as KL>0, K>N" + ?% and therefore L is always less than unity 
(as it must otherwise obviously be), and always positive. Clearly also, 
L increases with the increase of c and V, but diminishes with the increase of 
v, r, and N" (the last of which exists in K but cancels out from KL). 

We must now examine the conditions which hold respectively if / is to be 
of Type II or Type I, or is to be nearly or quite symmetrical. 

It will be of Type II, that is, always increasing in a drawn-out S-shape 
[(v) above] if L is not greater than -|- ; that is, if c$>v ~ Y + 2(N" -f r). But as 
just shown, in the cases which we are now considering, c must always be 
greater than v— -V-f-N-f r. Hence in order that / shall be of Type II, c must 
lie between the limits v— Y -f 2(N -f r) and v — Y + (N + r). Now if N + r is 
small, as is usually the case in short and sharp epidemics of zymotic diseases 
with long immunity, these limits will be narrow ; in other words, out of the 
whole range of possible values of c, that is of infectivity, only a small sector, 
and that the lowest one, will make/ of Type II. 

On the other hand, / will be of Type I, that is, of a bell-shape, regular or 
irregular, if L> | — that is, if c is any number greater than v—Y + 2(J8 + r). 
In this case /rises to a maximum ordinate \c and then falls to a constant 
level cL(l — L). Lastly, / becomes nearly symmetrical (see (v) above) when 
L = 1 nearly, that is, when c == v — V-f-(l + 7)(N" + r), where y is consider- 
able ;■ and also of course when 1ST and r are very small. These are probably 
just the conditions which hold in many of the short and sharp epidemics of 
zymotic diseases, such as measles, scarlatina, and dengue ; and they produce 
curves which differ from the perfectly symmetrical curve of c/K . dx/dt only 
by a small increasing term which, however, never exceeds c/K . (N-f-r)L. 
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For an example, suppose the case of a human infectious disease in which 
v—Y = 0*002, N" = 0*0004, r = 0*02, all assumed for a time-unit of one week. 
Then KL will not be positive unless c> 0*0224; that is, the happening 
will have no effect unless, say 1000 affected individuals infect or re-infect 
22*4 individuals per week on the average. If 1000 such persons can infect 
or reinfect 42*8 persons per week or less, the /-curve will be of Type II and 
will always ascend to the limit \c (that is, 0*0107 if c = 0*0428), at which 
figure it will remain indefinitely. If, however, c is greater than this — if the 
1000 affected persons can infect or re-infect say 100 persons per week — then 
the /-curve will rise to a maximum value of 0*025 of the total population, and 
afterwards decline to the constant value of 0*01647 of the total population, or 
16*47 per thousand, at which it will remain (L being 0*792 in this case). If 
c is still greater, if each affected person can affect one other person per week, 
so that c = 1, then L = 0*9795, that is, nearly the whole population 
ultimately becomes affected. Here /, the proportion of new cases, reaches a 
maximum of one quarter the whole population and thereafter declines to the 
constant value of 0*01997, or nearly 2 per cent, of the population, at which it 
remains indefinitely. (For further examples see Section IX, Part II.) 

(vii) Integral expressions for P and Z were obtained in equations 51 and 52 
by integrating both sides of equation 44. These may be written 

P = ? e vt (x/x . 6 ~ KL 'f- v)/K = Po(^7^o) (?; " v)/K 6^6-^- v ) L ^, (69) 

Z = Z 6 vt 6 KLt (x/x . €- KL ') c/K = Zo(x/x o y IK 6 vt 6~^'- Y ^ t . (70) 

Two important special cases present themselves here. If v— V = 0, the 
happening has no effect on the total number of the population, as in mild 
infectious diseases such as chicken pox or (?) dengue or party politics. 
This should be called the Equi variant Case of Proportional Happening. In 
it, P = P e w * as in Section IV (ii), and of course Z = Z e^. But if v—Y is 
small without being zero, care must be taken with the functions if c, N", and 
r are small. 

The other special case occurs when ^ = 0, and may be called the Case of 
Constant Natural Population. It is particularly important, since happenings 
which are considered during only a short period compared with the average 
length of life of the individuals concerned approximate to this type, because 
then v, the element of natural increase, is generally so small compared with 
the other constants that it may be neglected; in fact, the a priori curves 
of most epidemic diseases may be considered on this basis. 

If v = and V is negative (as in all injurious happenings) we see that, 

since 

P = P (^/^o)- v/(6 ' +V) 6 VLi , (71) 
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the population diminishes indefinitely as time progresses ; which is, of course^ 
what was to be expected since by supposition there is no natural increase 
to compensate for the loss due to the happening. 

If neither v = nor v— V — then, since x = L ultimately, that is, a 
constant, the question whether the population will finally increase or 
diminish depends upon the sign of the expression 

^__(^V)L = V + [(v-V)(N" + r)]/(c-v + V). (72) 

An important value of P is reached when the #-curve attains its centre 
of symmetry, that is when x = |L and t = r = 1/KL . log [(L— #o)/#o] 
(see (iii) above). Then 

P = P (L/2^ ) (r - v)/K [(L-~~^o)/^oF- (v - V)L]/KL ? (73) 

-Po(-|)^-^ K (LPo)^ L , 

if xq is very small, as when it equals 1/Po- 

Another important value of P is reached when Z = LP — 1, that is when 
i — 2r (equation 58). Then 

P = P [(L~:r )/^] [2 ^ ( ^ V)L]/KL . (74) 

A third important value of P is reached (if L>|) when the /-curve attains 
its maximum ; that is, when x = \ and t — r (equation 66). Then 

P = P (l/2^o) (v ~ V)/K [l/x .(L-x )/(2L- l)]C«-(«-v)L]/kl (75 ) 

These expressions are required to estimate the actual values of Z and F 
at these points ; but, if L is not very small, t will not be very large (see 
Section VIII (vii)), and therefore P can vary but little in that period from 
its original value P 0) and, at least for rough calculations, can be taken as 
remaining constant. 

(viii) It remains to consider the integrals of P, Z and F, which will be 
required for the analysis of several questions. "We proceed as in the case 
of constant happenings (Section IV (v)) by integrating both sides of the 
fundamental differential equations (44 and 45). Hence 

rt rt 

P--P == v Tdt-iv-Y) Zdt, (76) 

Jo Jo 

Z-Z = f W + (V-N-r)(W * (77) 

Jo Jo 



To interpret the first of these equations, we note that the actual change in 
the total population effected during the period t, namely P— 1\ equals what 

rt 

the change would have been by natural- variation only, namely vPdt, less 



j 



o 
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the difference between the natural- variation and the happening-variation in 

the affected part of the population, namely 1 (v-Y)Zdt. The second 

Jo 

equation, which is the same as equation 24, has been already interpreted in 
Section IV (v). It means that the total change in the number of the affected 
during the period t, namely, Z — Z 0) is equal to the sum of new cases, namely 

rt 

Felt, plus the affected immigrations, and less the affected emigrations, 

rt 
deaths, and recoveries, namely, (I— E — M— r)Zdt. Both these interpreta- 

jo 

tions agree with what was to be expected. 

The integrals of F and Z can be obtained at once in the important case 

when v = ; and are 

*t 
Zdt = (P-P )/V, (78) 

W = Z-Zo-(V-N"-r)(I > -P )/V. (79) 

Jo 

If in these we have also V — 0, then (P — P )/V becomes indeterminate since 
P remains constant; but the value of the fraction then is, from equation 71, 

( l» - P )/ V = Ii - i log, (L/x ). (80) 

(j 

The integrals of Z and F can also be easily obtained by direct integration 
of their values when v = 0. 

We should note that the total births, deaths, immigrations, emigrations, 

and recoveries among the affected during the period t are given, as with 

rt 
constant happenings, by multiplying Zdt by N, M, I, E, and r respectively. 



%> 



o 



The main result has already been summarised in the prefatory section. 
We can see from (vi) above that, for many values of the constants, the 
solution of the fundamental proposition stated in Section II yields curves 
generally similar to the curves frequently found in epidemics — that is, curves 
which are bell-shaped and nearly symmetrical, and tend to decline more 
slowly than they rose. This decline is in the case of the a priori curves due 
merely to the exhaustion of susceptible material ; but the further studies 
contained in the second part must be considered before attempts are made to 
apply the results to actual observations. 

For the projected contents of Part II see end of Section I above (p. 208). 
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